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by 

1/  2/ 

P.  C.  Tully-  and  J.  M.  Estes- 

ABSTRACT 

The  purpose  of  this  investigation  was  to  test  the  relative  consis¬ 
tency  of  the  vapor-liquid  equilibria  data  available  in  the  literature 
on  the  helium-nitrogen  system.  The  data  of  DeVaney,  Dalton,  and  Meeks; 
Kharakhorin;  Buzyna ,  Macriss,  and  Ellington;  Rodewald,  Davis,  and  Kurata; 
and  Gonikberg  and  Fastovskii  were  evaluated  by  an  isothermal  consistency 
test  equation  recently  proposed  by  Thompson  and  Edmister.  DeVaney1 s  data 
was  the  most  consistent;  Gonikberg' s  was  the  least  consistent. 

INTRODUCTION 

As  the  name  implies,  thermodynamic  consistency  tests  are  perform¬ 
ed  to  determine  whether  or  not  one  form  of  experimental  data  is  thermody¬ 
namically  consistent  with  another  form  of  experimental  data.  For 
isothermal  data  vapor-liquid  equilibria  data  are  usually  compared  with 
pressure-volume- temperature  data;  for  isobaric  data  vapor-liquid 
equilibria  data  are  usually  compared  WLt.h  enthalpy  data.  Since  there 
are  more  PVT  data  than  enthalpy  data  available  for  the  helium-nitrogen 
system,  the  isothermal  form  of  the  consistency  test  was  selected  as  a 
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Texas . 

2/  Research  Chemist,  Helium  Research  Center,  Bureau  of  Mines, 

Amarillo,  Texas. 


Work  on  manuscript  completed  October  1965. 


i 


■  j  ..  '  :  -  .•  .  •  .. 

'  :  . 

.  '  '  1  .  If  .  j  '  ,  .  •  ;  ■ 

•  '  -  •  ’  •  •  ■  ,  ■  •;;* 

.  -  j'  -  »:oj  - 

.  .  ' 

'*•••  '  {  .  ••  -  •  ... 

II  u 

•  •  U(j 

1  - •  .  :  /e 

i  •  ■ 

,  ’  i  i  ■  \  .  ' 

■  >  ...  -  ■  ■  ■  1 

.  •  t.  .  :  i  i  J  <  t-JV.'-. 

.  •  r, ••  .jr 


5 


basis  for  comparing  the  helium-nitrogen  vapor- liquid  equilibria  data 
available  in  the  literature. 

Consistency  tests  can  be  performed  in  either  the  differential 
form  or  the  integral  form.  The  former  is  a  point-by-point  evaluation 
of  the  data,  whereas  the  latter  averages  all  of  the  data  along  a 
given  isotherm.  The  integral  form  of  the  test  was  selected  because 
we  wanted  to  evaluate  the  data  of  each  investigator  in  as  few  segments 

as  possible. 

Five  sets  of  data  from  independent  investigators  were  selected 
for  performing  these  consistency  tests.  The  pressure  and  temperature 
ranges  of  these  data  are  given  in  table  1. 


TABLE  1.  -  Helium-nit.rogen  vapor- liquid  equilibria  data_eva_luat1ed 


Source. 

3/ 

DeVaney,  Dalton,  and  Meeks  (4)~ 
Kharakhorin  (7.) 

Buzyna ,  Macriss,  and  Ellington  (2) 
Rodewald,  Davis,  and  Kurata  (12.) 
Gonikberg  and  Fastovskii  (6) 


Pressure  range,  Temp,  range, 
psia  °  K  _ 

200-2000  77.0-120.0 

66-3152  68.0-111.5 

169-100.1  77.2-122.8 

200-1000  64.9-  77.2 

263-4196  78.0-109.0 


3/  Underlined  numbers  in  parentheses  refer  to  items  in  the  list  of 
references  at  the  end  of  this  report. 


The  investigations  of  Fedoritenko  and  Ruhemann  (5) ,  as  well  as 
those  of  several  other  authors,  have  not  been  included  in  this  group 
of  consistency  tests  because  they  did  not  publish  numerical  data. 

An  IBM  1620  digital  computer  was  used  extensively  in  carrying 
out  the  calculations  required  for  this  report. 
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In  K 


ISOTHERMAL  INTEGRAL  CONSISTENCY  TEST  EQUATION 
AND  ITS  EVALUATION 

The  consistency  test  equation  utilized  in  this  internal  report 
to  evaluate  the  helium-nitrogen  vapor- liquid  equilibria  data  is  one 
recently  proposed  by  Thompson  and  Edmister  ( 1^0  ?  which  is 


P 

1’  2 

x i  ( l  +  y151)d  ln  + 

p 

v i 


In  K? 

.  x2(1  +  yl5l)d  ln  K2 

ln  K2>P1 


P2  vL 


X1V1 


RT 


X2V2 


dP  (1) 


where 


(2) 


A  complete  derivation  of  equation  (1)  is  given  in  Appendix  I.  For  a 
vapor  which  obeys  the  Lewis  and  Randall  fugacity  rule,  § ^  —  0?  and 
equation  (1)  reduces  to  the  equation  of  Adler ,  et  al.  (jj)  . 

In  order  to  test  the  consistency  of  the  data,  it  is  necessary 
to  compare  the  summation  of  the  terms  on  the  left-hand  side  of  the 
equation  with  the  term  on  the  right-hand  side.  The  better  the  agree 
ment  is  between  the  opposing  sides  of  the  equation,  the  more  consistent 

the  data. 

The  quantities  x^,  y^,  K^,  x2 ,  y2 ,  K2 ,  T,  and  P  were  taken 
directly  from  the  data.  A  list  of  the  notation  used  in  this  report. 

is  appended . 

The  E  term  was  evaluated  using  an  equation  of  Thompson  which  was 

1 

derived  from  the  Bet  1.  in-  form  virial  equation  of  state.  His  derivation 


.  i.  nr  • 
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is  given  in  Appendix  XI.  While  the  original  derivation  includes  the 
second  and  third  virial  coefficients,  only  the  second  coefficients 
were  used  to  evaluate  in  this  report..  Thompson  ( 15)  says  this 
omission  does  not  significantly  affect  the  results.  The  virial  coef¬ 
ficients  used  were  calculated  using  an  equation  of  Wood,  et  ai„  (.16,)  • 

In  the  absence  of  complete  PVT  data  in  the  two-phase  region, 

the  volumetric  terms  on  the  right-hand  side  of  the  equation  were 

calculated  from  various  equations  of  state.  (f°r  the  mixture) 

was  assumed  to  be.  equal  to  (for  pure  nitrogen)  ,  and  the  latter 

V  -V 

calculated  by  an  equation  given  by  Strobridge  (.13).  an<^ 

were  calculated  by  the  Red lich-Kwong  ( 10)  equation  of  state,  and  a 
partial  molal  volume  equation  given  by  Prausnitz  (9) ,  respectively, 
was  calculated  from  the  partial  molal  volume  equation 

vv  =  yjl +  y2vi  (3) 

Use  of  these  equations  is  covered  in  detail  in  Appendix  III. 

Rodewald,  Davis,  and  Kurat.a  (.12)  published  values  for  the  den¬ 
sities  of  the  coexisting  phases  at  64.9,  69.3,  and  77.2°  K,  which 
were  obtained  by  using  a  system  of  simultaneous  equations.  These 
densities  were  inverted  to  obtain  the  molal  volume  terms  which  were 
used  to  evaluate,  their  data.  The  partial  molal  volume  term  was  ob¬ 
tained  from  the  Prausnitz  (9)  equation,  which  is  based  on  the  Red lich- 
Kwong  equation. 

The  integration  of  the  first  term  on  the  left-hand  side  of 
equation  (1)  was  carried  out.  by  evaluating  x^(l  +  y^)  and  in 
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at  each  pressure  along  a  given  isotherm  for  a  particular  investi¬ 
gation.  Least-squares  regressions  were  then  made  on  these  data 
using  the  x.  (1  +  y^)  term  as  ordinate,  and  the  In  K.  term  as 
abscissa,  and  second,  third,  fourth,  and  fifth  degree  polynomials 
as  models .  Each  of  the.  resulting  equations,  which  expressed 

x  (1  +  )  as  a  function  of  In  K. ,  was  integrated  between  the 

1  11  1 

limits  of  in  K.  evaluated  at  the  lowest  pressure  and  In  eval¬ 
uated  at  the  highest  pressure.  The  values  of  these  integrals  were 
used  to  select  the  polynomial  which  represented  the  data  best. 

When  the  integrals  of  two  consecutive  polynomials  agreed  to  within 
2  1/2  percent,  we  selected  the  lower  order  polynomial  for  use  in 
the  consistency  test. 

To  illustrate  this  procedure,  let  us  examine  the  polynomials 
of  the  first  integral,  A  ,  for  DeVaney's  77.0°  isotherm.  The  values 
for  the  first  integral  on  the  left-hand  side  are  tabulated  below. 


Degree  of  Polynomial 

A1 

Second 

.01637 

Third 

.01659 

Fourth 

.01657 

The  second  order  equation  was  selected  as  the  one  which  best,  repre 
sented  the  data  for  this  isotherm.  This  produced  an  equation  of 
the  form 

2 

=  A  +  B  In  Kt  +  C (In  Kj)  (A) 


xi(1  +  nfd 
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where  A,  B,  and  C  are  the  constants  determined  by  the  least  squares 
regression. 

The  second  term  on  the  left-hand  side  was  treated  in  much  the 
same  manner  as  the  first  term.  The  right-hand  side  of  the  equation, 
however,  could  not  be  fitted  with  an  equation  of  that  form.  There¬ 
fore,  an  equation  of  the  form 

VL  -  -  *2V2  =  A  +  ~P~  +  ~ 2 ~  ^ 

was  used. 

Typical  curves  for  the  three,  equations  are  shown  in  figure  2. 

The  integration  was  carried  out.  from  the  lowest  pressure,  point  1, 
to  the  highest  pressure,  point  2.  Some  investigators  added  the 
vapor-pressure  data  of  pure  nitrogen  to  their  vapor-liquid  equili¬ 
bria  data  to  test  the  consistency  of  their  data  using  other  consis¬ 
tency  test  equations.  This  technique  was  attempted  in  the  preparation 
of  this  report.  While  the  effect  on  the  terms  on  the  left-hand  side 
was  acceptable,  the  effect  of  adding  this  point  to  the  right-hand 
term  was  not.  From  figure  1,  it  can  be  seen  that  the  contribution 
made  by  the  pure,  nitrogen  point,  to  the  area  under  the  curve  is  very 
sizeable . 


'  b«  ...  t,  .  .  h,  a-j;.-,rfw 
_>jr-  i  i 

•  •  ...  -  .  ,  s>,  ' 

u  •  . 
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In  K, 


A  =  -.02490  = 


In  K  ,P2 

/  x^l+y^^d  In  Kl 


J  x2(l+y151)d  In  K2 


In  K^Pj. 


In  K2}Pl 


Figure  2. 


Typical  integrals 
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Furthermore,  the  shape  of  the  curve  between  the  two  points  has  a 
significant  effect  on  the  value  of  A3.  Therefore,  the  use  of  the 
pure  nitrogen  data  point  was  deleted  from  all  terms  to  eliminate 
the  introduction  of  this  uncertainty. 

DISCUSSION  OF  RESULTS 

The  results  of  the  consistency  test  are  given  in  table  2  and 
plotted  in  figure  3.  The  investigators  have  been  ranked  in  order 
of  their  consistency. 

Recalling  that  each  area  could  be  in  error  by  as  much  as 
2  1/2  percent,  the  most  probable  error  in  the  answer  would  be  4.3%. 
Therefore,  any  results  within  +  4.34  are  considered  consistent. 
Likewise,  if  any  two  investigators'  results  are  within  4 . 370  of 
one  another,  they  have  to  be  considered  as  being  of  the  same  order 
of  consistency.  This  is  illustrated  in  the  case  of  DeVaney  (4.6%) 
and  Kharakhorin  (7.0%). 

Examination  of  the  results  in  detail  indicates  that  the  high 
temperature  results  of  each  investigator  are  less  consistent  than 
the  low  temperature  results.  Virtually  none  of  the  data  above 
110°  K  is  consistent  at  all.  Four  possible  reasons  for  this  are; 

1.  The  Red lich-Kwong  equation  does  not  predict,  the  compressi¬ 
bility  of  nitrogen-rich  vapor  too  well. 

2.  The  liquid  sampling  methods  of  most  investigators  were 

such  that  fractionation  may  have  occurred. 

3.  The  values  for  may  be  more  in  error  at  the  higher 
nitrogen  concentrations  associated  with  the  higher  temperatures. 
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1/ 

table  2.  -  Results  of  consistency  te§£ 


Dif  f . 


T,°  K 

P  range, 
ps  ia 

Ai 

A2 

A1+a2 

A3  - 

A3  (a^+A2^  ( 

DeVaney , 
77.0 

,  et  a  1 . 
200-2000 

-.0164 

-1.65 

-1.67 

-1.54 

.  13 

80.0 

200-2000 

-.0196 

-1.66 

-1.68 

-1.52 

.  16 

85.0 

200-2000 

-.0249 

-1.47 

-1.49 

-1.48 

.01 

90.0 

200-2000 

-.0292 

-1.37 

-1.40 

-1.41 

-.01 

95.0 

200-2000 

-.0309 

-1.28 

-1.31 

-1.32 

-.009 

100.0 

200-2000 

-.0320 

-1.16 

-1.19 

-1.22 

-  .03 

105.0 

200-2000 

-.0300 

-1.00 

-1.03 

-1.09 

-  .06 

110.0 

400-2000 

-.0240 

-.482 

-.506 

-.554 

-.048 

115.0 

400-2000 

-.0142 

-.334 

-.348 

-.412 

-.064 

120.0 

600-2000 

-.0061 

-.100 

-.106 

-.155 

-  .049 

Kharakhor in 

68 . 0  66-*2 124 

-0.0122 

-2.59 

-2.60 

-2.46 

avr  . 

0. 14 

77.3 

71-2329 

-.0234 

-2.51 

-2.53 

-2.54 

-.01 

90.1 

71-3152 

- .0431 

-2.54 

-2.58 

-2.40 

.  18 

107.0 

514-2219 

- .0200 

-.338 

-.358 

-.465 

-.107 

111.5 

285-2866 

- .0178 

-.487 

-.505 

-.674 

-  .  169 

Buzyna , 
77.2 

et  al . 
169-1001 

-.0095 

-1.32 

-1.33 

-1.31 

avr  . 

-.02 

92.9 

169-998 

-.0183 

-1.04 

-1.06 

-1.17 

-  .  11 

112.1 

499-998 

-.0093 

-.190 

-.199 

-.239 

-.040 

118.1 

499-998 

-.0092 

-.158 

-.167 

-.153 

.014 

122.8 

499-998 

-.0106 

-.132 

-  .  143 

-  .082 

.061 

Rodewald ,  et  al. 
64.9  200-1000 

-  .0050 

-  .978 

-  .983 

-1.09 

avr  . 

-.11 

69.3 

200-1000 

- . 0066 

-.750 

-.757 

- 1 .  16 

-  .40 

77.2 

200-1000 

- .0093 

-1.20 

-1.21 

-1,15 

.06 

Gonikberg  and  Fastovskii 
78.0  263-4196  -.0326 

-1.22 

-1.25 

-1.43 

avr  . 

-.18 

90.1 

294-4053 

- .0384 

-  .949 

-  .987 

-1.19 

-  .20 

109.0 

391-4068 

.0051 

-.170 

-.165 

-.554 

-  .389 

7=  Diff. 
D^£— )(1007o) 


-8.5 

-10.3 

-.7 

.9 

.7 

2.1 

5.9 

8.6 

15.6 

31.9 

=  4.6  +  12.2 


-5.6 

0.5 

-7.9 

22.9 

25.1 


=  7.0  +  15.8 


-1.4 

9.6 

16.4 

-9.6 

-73.5 

=  10 . 7_ +  24.7 


10.8 
3^  .  6 
-5.0 

='U7T±  19.9 


12 . 5 
17.0 
70.2 


avr .  =  33.2  +  32.1 


1/  Rounded  from  standard  eight-digit  precision. 
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4.  Temperature  control  at  the  elevated  temperatures  may  not  be 
as  good  as  it  is  around  77°  K. 

Compressibility  factors  for  the  mixtures  were  calculated  using 
the  Red  1 ich-Kwong  equation  of  state  and  compared  with  the  experimen¬ 
tal  data  available  in  the  literature.  At  1,000  psia,  the  deviations 
were  4.1%,  3.1%,  and  3.4%  for  Rodewald's  64.9°,  69.3°,  and  77.2° 
isotherms,  respectively.  For  Canfield's  (3)  133°  isotherm  at  250 
atmospheres,  the  deviations  ranged  from  2%  at  y^  =  .15  to  7%  at 

y1  =  1-°- 

The  sensitivity  of  the  consistency  test  to  errors  in  x  and  y 
was  determined  for  DeVaney's  85.0°  and  120.0°  isotherms.  The  value 
of  x  was  reduced  by  107,  of  the  composition  or  0.01  mole  fraction, 
whichever  was  smaller.  The  value  of  y  was  increased  by  the  same 
amount.  For  example,  at  200  psia  and  85.0°  ,  x  was  changed  from 
0.0039  to  0.0035;  the  value  of  y  was  increased  from  0.806  to  0.816. 
These  changes  produced  an  increase  in  the  K-value  of  12.87,  (206.7 
vs.  233.1).  Similar  perturbations  were  made  at  all  points  along 
the  two  isotherms.  For  the  85.0°  isotherm,  the  consistency  percent 
difference  changed  from  -.7  to  -18.7.  For  the  120  isotherm,  the 
percent  difference  changed  from  31.9  to  27.2. 

The  magnitudes  of  the  three  volume  terms  which  go  to  make  up 
the  right-hand  side  of  the  consistency  equation  were  also  investiga¬ 
ted.  Values  of  VL,  ,  and  are  given  in  table  3  for  DeVaney's 
77.0°  and  120.0°  isotherms  at  the  lowest  and  highest  pressure  run 
on  each  isotherm.  Examination  of  the  penultimate  column  indicates 


ni 

. 


TABLE  3 .  -  Magnitude  of  terms  on  right  side 


T, 

°  K 

P, 

ps  ia 

Xi 

77 

200 

.0026 

77 

2000 

.0214 

120 

600 

.0253 

120 

2000 

.140 

>  -» 

l> 

x2 

vv 

2 

481 

.997 

419 

62.2 

.979 

54.0 

427 

.975 

97.5 

98.9 

.  860 

46.0 

for  DeVanev's  77°  and  120°  isotherms 


VL  -(X^CV/ 


(Vl)V 

-(x^^ 


(V'-(x1)(v,)v 

-  (Xj)  (V2)V) 
RT 


34.6 


-1.25 


-418 


34 . 6 

-1.33 

-52.8 

53.1 

-10.8 

-95.0 

53.1 

-13.8 

-39.6 

-384  -.00414 

-19.5  -.000210 

-52.7  -.000364 

-.3  -.00000021 


•• 

V>  \  — i 
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that  the  value  of  the  right-hand  side  of  the  equation  is  largely 
attributable  to  the  low  pressure  data.  This  minimized  the  effect 
of  the  error  induced  by  assuming  that  the  liquid  is  pure  nitrogen. 
The  magnitude  of  is  seen  to  be  important  at  all  points.  Thus, 
the  consistency  test  is  dependent  on  the  ability  of  the  Redlich- 
Kwong  equation  to  predict  the  compressibility  of  the  vapor. 


CONSISTENCY  TEST  USED  BY  BUZYNA  AND  RODEWALD 

Both  Buzyna  (2)  and  Rodewald  (12)  used  Adler  s  consistency 

test  on  their  data.  Neither  of  these  investigators  incorporated 

the  £  term  to  account  for  the  nonideality  of  the  vapor.  Both 

1 

used  graphical  integration  which  is  less  accurate  than  the  digital 
computer  methods  utilized  in  this  report.  Both  used  the  vapor  pres¬ 
sure  of  nitrogen  as  the  lower  limit  of  integration,  which  we  did  not 
Both  of  these  investigators  made  simplifying  assumptions  to 
carry  out  their  consistency  tests.  Buzyna  made  the  assumption  that 
=  V^.  He  used  four  combinations  of  sources  of  volumetric  data 
and  got  widely  differing  answers  by  these  four  methods.  Two  of  the 
sources  involved  the  use  of  generalized  charts.  The  other  two  were 
Otto's  (8)  equation  for  helium  compressibility  factors  and  experimen 

tal  data  for  Z^.  Even  when  he  used  Otto's  equation  for  Z. ,  he  still 

Z  1 

used  generalized  charts  to  compute  Z  . 

Rodewald  used  his  own  volumetric  data,  of  course.  He  assumed 

that  zy  =  ZV.  Thus,  he  was  able  to  simplify 
1  m 

zv.  .  ,  ... 

At  the  low  temperatures 


T  L  V  /  1  1  \  ml 

LZm  +  Z1  Y  1\K  "  K.)  “  kJ 


r  l 

vi 

to 

z 

-  z 

_  m 

m_ 
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where  his  data  were  taken  (64°  to  77°),  this  is  not  as  serious  an 
error  as  it  is  at  120°  (y2  increases  with  temperature)  . 

There  is  qualitative  agreement  between  the  consistency  test 
results  of  these  two  investigators,  and  the  results  of  this  report. 

CONCLUSIONS 

Before  the  results  of  this  work  can  be  made  conclusive,  it 
will  be  necessary  to  correct  some  of  the  deficiencies  indicated 
above  in  the  prediction  of  thermodynamic  properties.  Other  consis¬ 
tency  tests  can  also  be  performed.  This  work  is  currently  underway. 


.. 

,  t ::  .. 

■  ■  ■  '  ti  :  ' 
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APPENDIX  I 

DERIVATION  OF  THOMPSON -EDMISTER  EQUATION 
This  test  is  based  on  the  following  general  form  of  the  Gibbs- 
Duhem  equation 

S  n .  d  (j, .  =  -SdT  +  VdP 

1  1 

This  equation  can  be  derived  as  follows. 

Consider  first  a  reversible  process  acting  on  a  closed  system 
wherein  fluid  pressure  is  the  only  acting  force,  then  by  the  First 
and  Second  Laws  of  Thermodynamics 

H  =  U  +  PV  (I-D 

dH  =  dU  +  PdV  +  VdP  (1-2) 

U  =  TdS  -  PdV  (1-3) 

Substituting  the  value  of  dU  from  equation  (1-3)  into  equation 
(1-2)  gives 

dH  =  TdS  +  VdP  (1-4) 

The  Gibbs  free  energy  function  G  is  defined  by 

G  =  H  -  TS  (1-5) 

from  which 


dG 


dH  -  TdS  -  SdT 


(1-6) 


>• 


* 


22 


Substituting  the  value  of  dH  from  equation  (1-4)  into  equation  (1-6) 

d G  =  VdP  -  SdT  (1-7) 


In  order  to  introduce  where 


l  P,T 

now  consider  an  open  (variable  mass)  system,  where 


(1-8) 


G  =  G(T,  P,  nl,  ...  nn) 


(1-9) 


General  differentiation  gives 


From  equation  (1-7) 


(1-11) 


-S 


(1-12) 


Substituting  these  values  in  equation  (1-1.0) 


dG 


-SdT  +  VdP  + 


I  (Sr) 


i  T,P 


dn. 

i 


(1-13) 


Specifying  conditions  of  constant  P  and  T,  and  integrating  gives 


G  =  E  u, .  n  . 

i  i 

(See  equation  (1-8)  for  jj,  .) 


(1-14) 


'■ 
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General  differentiation  of  equation  (1-14)  gives 

dG  =  Hu.  dn.  +  E  n.  du .  (1-15) 

ii  i  i 

(At  conditions  of  constant  P  and  T,  E  n_^  dp,^  =  0.) 

Comparing  equations  (1-13)  and  (1-15) 

En.  da.  =  -SdT  +  VdP  (1-16) 

l  r  l 

This  completes  the  derivation  of  the  Gibbs-Duhem  equation. 

Equation  (1-16)  will  now  be  written  for  a  closed  system  in 
equilibrium.  Dividing  equation  (1-16)  by  the  total  number  of  moles, 
E  n.  where 

l 

x  .  =  n . /E  n . 

l  ii 


E  x.  du.  =  -SdT  +  VdP 

i  'i  —  — 


(1-17) 


where  is  moial  entropy  and  V  is  molal  volume.  Since 
(U^  —  iu.(P,T,x.) 


From  equation  (1-11) 


d_ 
LdP 


dG 

dn 


P,T  T,x 


,x.  P  ,T 

i 


P,T 


V. 

l 


(1-35) 


>• 


' 


. 
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Similarly,  from  equation  (1-12) 


Substituting  in  equation  (1-18)  gives 


V.  dP 

l 


S.  dT  + 

l 


T,P 


(1-19) 


Combining  equations  (1-17)  and  (1-19)  at  constant  P  and  T, 


x  •  \  > 

1  ^Sxiyp,T  1 


dx .  =  0 


(.1-20) 


Equations  (1-19)  and  (1-20)  can  be  written  separately  for  liquid 
and  gas  phases,  or 


.  L 

du, .  = 

1 


dq,V 


VL  dP 

i 


-  SL  dT  +  (~^)  dx. 

1  \OX./B  rn  1 

1  P,T 


5pV 

VV  -  SV  dT  +  1  1 

i  i 


By7P,T  dyi 


(I-21a) 


(1.-2  Lb) 


I-  (r1 


i  P  ,T 


dx .  =0 

i 


(I -22a) 


Yy.  A  dy.  -  0 

A  iAsVp,t  1 


(I-22.b) 


If  the  phases  coexist,  in  equilibrium,  an  additional  condition 
is  imposed  which  is 


-  < 


(1-23) 


' 
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Substituting  equations  (I-21a)  and  (I-21b)  into  equation  (I  23) 


and  rearranging 


(V  -  V^)dP  -  -  S^dT  (1-24) 

AV  dP  -  AS.  dT  (1-25) 

i  i 


where  AS  =  partial  molal  entropy  difference  between  equilibrium 
i 

vapor  and  liquid. 

Every  term  in  equation  (1-25)  is  multiplied  by  x..  The  resulting 
equation  is  expanded  for  a  binary  system 


+  x2AV2 


dP 


+  x2^2 


(1-26) 


so  that  equation  (1-26)  describes  the  liquid  phase  of  a  binary  system 
Expanding  equation  (I-22a)  for  a  liquid  phase, 


Bp,. 


l\Bx 


1  P  a  T 


dx  X2 


Bp, 


Bx 


dx. 


=  0 


2  P  jT 


(1-27) 


and  (I-22b)  for  the  vapor  phase 


0 


(.1-28) 
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Rearranging  and  recalling  that  dy^  =  -dy^ 


(S 


2  P  ,T 


a  V 

II  (itT) 

y 2  '^y y p  'p 


(1-29) 


Substituting  equations  (1-27)  and  (1-29)  into  equation  (1-26)  and 
collecting  terms 


V 


.y -x.,  ,d|ilN  _  _ 

(iV  )(^)  dyl  =  (X1A?1  +  W  dP  -  (xlAV  W  dT  (I-30) 


^1  7  p  T 


dp, 


V 


/u^]\ 

An  expression  will  now  be  derived  for  -  BY  definition  of 


1  P,T 


the  fugacity  f  at  constant  T 

dG  =  RT  d  In  f  (dT  =  0) 


(1-31) 


Therefore 


duV  =  RT  d  In  fV  (dT  =  0) 

pi  i 


(1-32) 


and 


V  -V  V 

u.  =  RT  In  f .  +  p  . 
*1  l  i 


(1-33) 


from  which 


(S1). 


,d  In  fV 


(1-34) 


From  equation  (1-11) 


.*Y. 


\ 

r 9 

_ 

l> 

H 

> 

so 

T  ,n . 

Ldn. 

l 

T ,n . 

T,P 

'^ni  i  ,P 

(1-35) 
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Substituting  equation  (1-35)  into  equation  (1-34) 


-V 
In  f\ 

hif-1). 


Integrating 


In 


(T  constant) 


If  state  1  is  chosen  such  that  f,  -  Py^,  then 


In 


where  a. ^ 


(1-36) 


(1-37) 


(1-38) 


(1-39) 


Expanding  and  rearranging  equation  (1-36) 


RT  In  f.  -  RT  In  P,  -  RT  In  y. 

i2  1  i 


RT  In  P2  -  RT  In 


(1-40) 


P 


1 


or 


RT  In 


(1-41) 


Letting  P^  -♦  0  and  P2  -»  P 


. 


-• 

. 
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RT  In 


fV 

1 

Py, 


o 

0 


a.dP  = 

l 


-  j  (?  -  '> 


0 


from  which 


RT  In  fV  =  RT  In  y.  + 

l  i 


V.dP 

l 


0 


Substituting  equation  (1-43)  into  equation  (1-33) 


=  RT  In  y.  +  J  V.dP  +  n.V 


0 


Differentiating 


§■) 


p  >av4. 


i  P jT 


RT 

y a  *<  i#  p  t 

1  0  ’ 


dP 


If  we  let 


X. 

RT 


P  — 

Lav, 


dP 

o'dy./p  p 

0  5 


then 


'•  (r-1)  =  RT  (1  +  y .5  .) 

1  VaVp,T 


(1-42) 


(1-43) 


(1-44) 


(1-45) 


(1-46) 


(1-47) 


The  equilibrium  distribution  ratios  (K-values)  will  now  be 
incorporated .  From  the  definition  of  K-value 


v 


— 


. 


•  . 

b.  t 
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K 


1  x . 

1 


(1-48) 


it  follows  that 


In  =  In  -  In  x^ 


(1-49) 


Differentiating  and  multiplying  both  sides  by  x 


X1  d  yl 

x.  d  In  K.  -  -  dx. 

1  1  y,  1 


(1-50) 


Similarly  for  the  other  component 


x2  d  y2 

x_  d  In  K  =  -  dx, 

2  2  y0 


(1-51) 


Adding  equations  (1-50)  and  (1-51) 


x  d  y  x2  d  y2 

X1  d  In  Kx  +  x2  d  In  K2  =  —  +  “ 


(1-52) 


-  (- 


yl  ~  xlVdyl 

i  -  y,Ayi 


(1-53) 


The  chemical  potential  and  K-value  expressions  will  now  be 
incorporated  into  the  equilibrium  equation.  Substituting  equations 
(1-53)  and  (1-47)  into  equation  (1-30) 


x^  d  In  +  x2  d  In  K_2 


(x  +  x2AS2)dT  -  (x][AV1  +  x2AV2)dP 


RT(1  +  y^p 


(1-54) 


i  i  .  .  i  .  }  h 


r  ,<  .  o 
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Restricting  equation  (1-54)  to  an  isothermal  path, 


S  In  K  ^  -2\ 

\rsr~L +  x2\  sp  ) 


b  In  K, 


x^AV^  +  x^/XV^ 

RT(1  + 


(P-55) 


VL  -  -  x/2 

RT(1  +  y^p 


(1-56) 


can 


Every  term  can  be  multiplied  by  the  pressure  and  the  equation 
then  be  expressed  in  compressibility  factors  rather  than  volumes 


PV  yV 
RT  5  1 


PVL  -v  .  fA 

RT  ,  and  Z2  RT 


Recalling  that 


dP  d  In  P 


x . 


b  In  K. 
l\d  In  P 


+  x. 


B  In  K, 


2\d  In  P  /  j. 


ZL  -  x  ^  -  x2 
(l  +  y^j) 


(1-57) 


The 


numerator  of  the  right  side  can  be  expressed  in  values  of  y  ^ 


K  ,  and  rather  than  liquid  compositions 

12 


ZL  -  X]2[ 


*2*2 


=  Z 


2\ 


y. 

v  ^ 


(1-58) 


V  -V 

Z  Z1  yl  -V 

Substituting  ~  ~  for 

y  o 


ZL  -  : 


,  <  -  ^  ^  ■  @  ^  -  ex 


V  -V 

yoX/ZV  -  Zl  yl 


y , 


(1-59) 


vrV  f  1  1 


“  Z  '  yl  Z1  \K.  K, 


Z 

K 


V 


(1-60) 


' 


1  ■  "  <■  -‘I 
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Using  this  expression,  equation  (1-57)  becomes 

/S  ln  K1\  ^  ln  K2^  ZM1x  ~  y],Zl(1-/Kl  ~  1/K2)  '  ZMix/K2 

XiU  In  p'7t  +  X2\9  In  P  =  U  +  y^) 

Equation  (1-61)  can  also  be  expressed  in  integral  form,  as  used 
by  Thompson  and  Edmister  (14) 


(1-61) 


In  ,  P2 

I  xi(1  +  ypi)d  ln  Ki 

ln  K  ,  ?l 


In  K2 ,  P2 


+ 


x2(l  +  y1§1)d  in  K2 


In  K2,  Pl 


■  JD 


JMix 


+  y- 


Z^(1/K2 


1/K.p 


zL.  /K 

Mix  2J 


d  in  P  (1-62) 


If  the  Lewis  and  Randall  rule  holds  for  the  vapor,  ^  =  0 
and  equation  (1-62)  reduces  to  that  derived  by  Adler,  et  al. 
Equation  1-56  can  be  expressed  in  integral  form 


In  K  ,  P2 

Xf(i  +  y1?1)d  ln  Ki  + 

in  Kl,  Pl 


In  K2 ,  P2 


x2(l  +  y1?1)d  In  K2 
ln  K2 ,  P^ 


P 


1 


-V  -V 

XjVj  -  x2V2 

RT 


(1-63) 


Equation  1-63  is  the  version  actually  used  in  this  consistency  test. 


' 

■ 
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APPENDIX  II 


VIRIAL  EXPRESSION  FOR 

Thompson  (15)  has  derived  an  expression  for  ^  using  the  Berlin 
form  virial  equation.  Following  is  Thompson's  derivation 

(11-1) 


V 


Y  +  B7  +  c'p 


and 


V  =  (n^+n^) V  - 


(n.+n„)RT  { 

- - - - -  +  (nL+n2)B  +  (n1+n2)C  P  (H-2) 


The  Berlin  form  second  and  third  virial  coefficients  are  expressed  in 
Leyden  form  equivalents  because  of  their  preponderance  m  the  litera¬ 


ture 


B 


B  =  yiBn  +  2yly2B12  +  y2B22 


(II-3) 


n, 


2n  n„ 

B,,  +  - B .  „  + 


n. 


B, 


(n^)2  U  (i^+np2  12  (n^)2  22 


(II-4) 


C-B 

RT 


L 

RT 


n. 


(nl+n2) 


3  Clll  + 


9  2 
3nin2 

(n  +n2) 


3  112 


+ 


9  2 
3n^n2 

(ni+n2) 


3  122 


n. 


+ 


(n1+n2) 


3  222  J 


RT 


n. 


4 


^nl+n2^ 


2  2 

2  4nln2  2 

-  B7.  +  - L-L—-  B“ 

4  il  ,  ,  N4  12 

(n1+n2) 


n, 


4  /  3 

2  4nln2 

b:„  +  — ■ — 


2  2  ,  3 

2n  n  /  ^n]_n2 

B ,  ,B1fl  + - -  ~  r  B,„B_  + 


+  4  b22  +  x4  d11u12  t  ,  ,  .4  "11"22 

(n,+n,)  (n,+n,)  Oj+V 


1  2' 


1  2' 


(Vn2> 


<-B12B22')  J 


(II-5) 
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Substituting  equations  (II-4)  and  (II-5)  into  equation  (II-2) 


(n  +n  )RT  p  n2  2n  n  n? 

V  =  P  +  Ln1+n2  B11  +  nL+n2  B12  +  B22j 


P  r  ni 
+  -L— 


2  2 
3n  n  3n  n 

c, , .  +  — -  ^  c,  ,„  +  — — —  c 


RT  L,  ,  ,2  111  '  ,  .2  ~112  '  ,  s2  ~122 

nl+n2^  (n1+n2)  (n  +n2) 


n. 


+ 


2  2  4 

_2  4nln2  2  n2 

B,  ,  +  - r  B.  „  +  - — 


B, 


(  ,  2  222J  RT  L  ,3  "11  '  ^  3  “12  '  ^  ,3  ~22 

^  1+  2>  (n1+n2)  (n]L+n2)  (np-np 


4nln2 

+  - —  —  B ,  ,  B ,  „  + 


2  2  3 

2nn  4nn 

- ^ —  B,  ,B^  +  - B,  B 


,  v3  11  12  '  ,  .3  ~11~22  '  ,3  ~12~22j 

(np-np  (n1+n2)  (n  +n) 


(II-6) 


The  definition  of  §  ,  equation  2,  contains  the  partial  molal  volume 

2 


V. 


/ 3V  \ 

RT  ,  r 

\dn  /  ^ 

"  p  r  _ 

1  P,T,n 


2  2 
r_nl+2nin?  2n? 

- 2  Bn  + - 


n. 


(ni+n2) 


(np-n) 


2  B12 


(nl+n2) 


2  B22 


+ 


(ni+3nin2) 


6n  n2 
+  - C 


RT  L,  ,  >.  3  wlll  '  ,  ,  v  3  ~112 

(nl+n2)  (ni+n2) 


(3n2-3n1n2) 

(nl+np  3 


122 


2n; 


(ni+n2) 


3n, 


(nx+n2) 


3  122 J  RT 


r(n4+4rutO 


l'’“l“2/„2  .  (8nln2  '  4nln2)  2 

4  BU  +  ,  .  .4  B12 


(nl+n2) 


(ni+n2) 


19  2  2  /  3  .  2  2 

2  12n  n  4n  n  -2n  n 

J  B  +  - B  B  +  B  B 

4  22  ,  .  v4  11  12  ,  .  v4  ii  22 


(nl+n2) 


(ni+n2) 


4n2-8n  v? 

+  — - - ^  B  B. 


(n1+n2) 


4  12  22 J 


(I I- 7) 


(  !-  I  •  ■.!  !  "  *■  i.  .  it?. 


. 
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Dividing  by  RT,  subtracting  1/P,  and  expressing  the  coefficients  in 
mol  fractions 


Ii  _l 

RT  P 


RC  [(yl+2yly2)Bll  +  2y2B12  -  y2B2 2]  +  72~2  [(yl+3yly2)Clli 

K  1 


+  6yly2C112  +  ('3y2"3yly2^>C122  '  2y2C222_  ~  2  2  Lyl+4yly2')Bll 

K  1 


+  (Sy^-^iy^L  ■  3y2B22  +  12yly2BllB12  +  (4yiy2“2yiy2} 


(‘BilB12')  +  (‘4y2"8yly2')B12B22J 


(II-8) 


Integrating  and  expressing  concentration  in  exclusively 


(rt  -  ~f)‘ 


vRT  P7dP  RT 


L(2y1-y2)Bu  +  U^y^y^B^  -  (l^y.+y^B^ 


+ 


2R  T 


2  2  L-2r2y3i)cm  +  (6yi-12yi+6y?)cii2  +  (3-12yi  +  15yi 


-  6yi>C122  -  (2-6yi+6y2-2y3)C222]  - 


2R  T 


oy3-3y4)B2u 


+  (8yi-28y2+33y3-12y^)B22  -  (3-12yi+18y2-12y3+3y4)B22 


+  (12y2-24y3+12y^)BllBl2  +  Oy^+iey^yX^ 


+  (4-24y1+48y2-40y3+12y4)Bl2B22 


(II-9) 


Equation  2,  the  definition  of  can  be  written 


P 


S,  = 


dy 


1 

0 


/V  -  x 

U  -  -±r) dP' 

\RT  P  J 


(II- 10) 


1 


. 
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Performing  tine  indicated  differentiation  and  expressing  concentration 

in  y  where  convenient 
7  2 


^1  RT 


_2y2Bn-4y2B12  +  2y2B22_ 


+ 


2  2 
2R  T 


6yiy2Clil+6y2(y2“2yi)Cil2 


“  6y2(-2y2  yi^C122  +  6y2C222. 


2 

2R  T 


i125'ly2Bll+8y2(1‘6yly2)B12 


+  12y2B22  +  24y1y2(y2-y1)BuB12  +  4y2(i-6y1y2)B.L1B22 

-  24y2(y2-y1)B12B22^ 


2Py, 


r  q 

3P  y2  r 

B  +B„„-2B , _ 

L  11  22  12_ 

2  2 

R  T 

LylClU+(y2'2yl)GU2  +  (yl'2y2)C122 


+  y0c 


2 

2P  y 


2  2„2  .  „  2^2  -  .  — 2 


2  222J  r2t2 


3y1Bu+3y2B22+(l-6y1y2)(BuB22+2B12) 


+  6(y2'yl)<:ylBllB12'y2B12B22). 


(11-12) 


2Py, 

~rt" 


B. ,+B,,  -2B 


11  22  12 J  2RT 


+ 


3P 


[yicm+(y2-2yi)c 


112 


+  (y1-2y2)c122  +  y2c222_ 


RT 


'  t  2  J2  2  2  . 
3  'ylBll+y2B22^ 


+  (l-6y1y2)(BuB22+2B22)  +  SB^-y^  (y1B11-y2B22>]) 


(11-13) 


Attempts  to  derive  an  expression  for  from  the  Leiden  form 

virial  equation  of  state  were  unsuccessful,  as  were  attempts  to 
derive  an  expression  from  the  Rediich-Kwong  equation. 
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APPENDIX  III 

The  Redlich-Kwong  (10)  equation  of  state  is 


P 


RT 

V-b 


_ a  _ 

T°“5  V (V+b) 


(III- I) 


It  also  takes  the  form 


1  A2/  h  \ 

1-h  B~\l+h/ 


(III-2) 


where 


2  2.5 

R  T 


?  5  2  5 

0.4278  T  P  r  ’  3 
c  c 


(III-3) 


0.0867  T 

_b_  =  c 

RT  P  T 

c 


(III-4) 


(III-5) 


Equation  (III-2)  was  programmed  for  the  computer.  This  program 
solves  for  Z  by  an  iteration  scheme  involving  interval  halving. 
This  program  used  the  mixing  rules  suggested  by  Redlich  and  Kwong , 


i .  e . 


(Ill- 6) 


B 


(III- 7) 


J 


f 
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2  2 

3  =  ai  yl  +  a2  y2  +  * 


•  +  2a 12  yl  y2  + 


(III-8) 


a . 


12 


(a L  a2) 


1/2 


(III-9) 


A  =  >  y.  A. 

j  j 

j 


(III-10) 


The 


:V  . 


Prausnitz  ( 9  )  equation  for  is  based  on  the  relationship 


V. 


-(dP/dn  ) 

1  1 5V 3n^ 

(SP/dV)„ 

/T,n1,n2 


(III- 11) 


The  corresponding  equation,  which  he  derived  from  the  Rediich-Kwong  is 


vi  ■ 


T1/2V(V+b) 


{2yi  al  +  2y2^ai  a2^"  La  b1/^V+b)  }  ‘ 


RT/  (V-b) 


1  + 


b  /(V-b) 

•  J- 


(a/T1/2) 


(2V+b)/V2  (V+b)2 


RT/ (V-b) 


(III- 12) 


The  Strobridge  (13)  equation  for  the  density  of  saturated  liquid 


nitrogen  is 


p  =  k  +  k  x  +  kQ  x2  +  k.  x  +  k  x  (III- 13) 

H  1  2  3  4  5 


where 


x 


126.26 


1/3 


(III  - 14) 


The  constants  k^  through  k^  have  numerical  values  given  by  Strobridge 


■ 


' 
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APPENDIX  IV 


At  the  beginning  of  this  test,  most  isotherms  had  equations  as 
high  as  the  fifth  degree  fitted  to  them  and  some  went  as  high  as 
the  thirteenth  degree.  These  higher  degree  equations  frequently 
gave  better  apparent  fits  to  the  curve  (i.e.,  smaller  residuals). 
However,  integration  of  these  equations  produced  impossible  results. 

For  the  right  side  of  the  equation  for  one  isotherm  twelve 
different  polynomial  curves  were  fitted.  They  were  of  the  type 
VL  -  XjV,  -  x/2  =  A  +  ~~  +  -~- 


=  A'  +  ^  +  -c:  +  -Bi 

P  2  3 

p  p 


etc  .  ,  up  to 


=  A"  H — H — -I — r- 


M 


.13 


P  P  P 

Table  IV- 1  lists  the  degree  of  each  of  these  polynomials,  the 


standard  estimate  of  error  for  each  one  (based  on  how  well  the  curve 


fits  the  data) ,  and  the  area  under  each  of  the  integrated  polynomials 


TABLE  IV -1  -  Integration  of  various  polynomials 


>f  Polynomial 

Standard  Estimate 
of  Error 

Area  Under 

Integrated  Polynomial 

2 

.000702 

-4.00 

3 

.000474 

-4.02 

4 

.000269 

-3.68 

5 

.0000875 

-11.05 

6 

.0000922 

6.77 

7 

.0000164 

-15,844 

8 

.0000178 

187  ,360 

9 

.0000197 

-1,231,247 

10 

.0000219 

1,252,326 

11 

.0000252 

1,466,312 

12 

.0000311 

1,686,911 

13 

.0000384 

7,247,713 

* 
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The  first  two  answers  were  in  good  agreement.  The  next  was 
poor,  and  the  rest  absurd.  It  developed  that  the  higher  degree 
equations  were  producing  "tear  drops"  between  the  actual  points. 
This  is  illustrated  in  figure  IV- 1. 


T 

>  CM 

l> 

CM 

X 


>  1 
l> 

x 

i 

h-q 

>1 


A  + 


actual  points 
reasonable  curve 
curve  fit  from 

B  ,  C_  1_ 

2  •  •  •  • 


p 


+ 


8 


Figure  IV- 1.  -  "Tear  drop"  from  higher  degree  equation  for  A, 


The  two  lower  degree  equations  produced  no  such  "tear  drops." 

Thereafter,  only  second,  third,  and  fourth  degree  equations 

were  fitted  for  A^  and  A^  ,  For  A^  equations  of  the  first.,  second, 
and  third  degree  in  the  reciprocal  of  P  were  generated. 


>' 


it  >.  > 


\ 
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NOMENCLATURE 

A 


A  = 

*  — 

=  parameter  of  Red lich-Kwong  equation 

=  area  of  an  integral 

a  = 

:  parameter  of  Redlich-Kwong  equation 

B  = 

=  parameter  of  Redlich-Kwong  equation 

=  second  virial  coefficient 

b  = 

=  parameter  of  Redlich-Kwong  equation 

C  = 

=  third  virial  coefficient 

f  = 

=  fugacity 

G  = 

=  Gibbs  free  energy 

H  = 

=  enthalpy 

h  = 

=  parameter  of  Redlich-Kwong 

BP 

equation  =  ~ 

K  = 

-  vapor-liquid  equilibria  phase 

distribution  ratio,  y/x 

=  Kelvin  temperature 

k  = 

=  constant  in  Strobridge  equation 

n  = 

=  number  of  moles 

=  number  of  components 

"  P 

=  pressure 

R 

=  universal  gas  constant 

S 

=  entropy 

T  = 

=  temperature 

U  = 

=  internal  energy 

V  = 

=  volume 

x  =  liquid-phase  composition 

y  =  vapor-phase  composition  in  equilibrium  with 
Z  =  compressibility  factor,  PV/RT 

Greek  Letters 

O'  =  RT/P  -  V 

A  =  change  in  a  quantity 

(j,  =  chemical  potential 

p  =  density  in  gram  moles/liter 


Subscr ipt s 

c  =  critical  property 
.  th 

1=1  component 
.  th 

j  =  j  component 
Mix  =  mixture 

1  =  lighter  component 

=  first  integral  of  consistency  test 
=  lower  limit  of  integration 

2  =  heavier  component 

=  second  integral  of  consistency  test 
=  lower  limit  of  integration 

3  =  third  integral  of  consistency  test 

1-5  =  different  constants  in  Strobridge  equation 


subbar 


molal  quantity 


I  -  '  :  lU»lO 


: 


Superscripts 


L  =  liquid  state 
V  =  vapor  state 
°  =  ideal  state 

=  Berlin  form 

superbar  =  partial  molal  quantity 

Abbreviations 

In  =  logarithm  to  the  base  e 
PVT  =  pressure-volume- tempera ture 


R-K 


Red lich-Kwong 


r.  { 

■ 


A 


i1 

f 


